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Abstract. Cir X-l is one of the most intriguing galactic X-ray sources. It is a ~ 16.6 days variable X/radio source, 
a type-I X-ray burster and a QPO emitter. In spite of an uncertain optical counterpart classification, all these 
properties identify the source as an LMXB. The morphology of the surrounding radio nebula has suggested an 
association with the nearby (~ 25 arcmin) SNR G321.9-0.3, implying that Cir X-l is a runaway binary originated 
from the supernova explosion ~ 10 years ago. To investigate this hypothesis, we carried out a proper motion 
measurement of the Cir X-l m ~ 19 optical counterpart using a set of HST/WFC and WFPC2 observations taken 
~ 8.6 years apart. We obtained a 3<r upper limit of w 5 mas yr _1 on the source proper motion. Since the runaway 
hypothesis would have implied a proper motion due North ranging between 15 and 75 mas yr -1 , depending on 
the actual age of the SNR, our result definitively rules out the association between Cir X-l and SNR G321.9-0.3. 
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1. Introduction. 

Cir X-l is a strong and persistently bright galactic X- 
ray source with a flux modulated over a 16.6 days cycle 
(Kaluzienski et al. 1976, Bignami et al. 1977), which is 
commonly believed to coincide with the orbital period of 
a binary system in a highly eccentric orbit (Murdin et 
al. 1980). Type-I X-ray bursts have been discovered by 
EXOSAT (Tennant, Fabian and Shafer 1986a,b), as well 
as high/low frequency QPO (Tennant 1987,1988) whose 
properties suggested that Cir X-l might represent an ideal 
link between the classes of the canonical Z and atoll-type 
sources (Oosterbroek et al. 1995). In addition, long-term 
variations in the X-ray flux have been observed, probably 
due to a partial obscuration of the central source by the 
precession of a warped accretion disk (Oosterbroek et al. 
1995). Recently, X-ray spectral evolution of Cir X-l along 
the orbital cycle has been reported by Iaria et al. (2001). 
While the general X-ray properties of Cir X-l (e.g. type-I 
bursts, QPOs) point towards a LMXB, its nature is still 
unclear. Although an optical counterpart has been identi- 
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fied (Moneti 1992; Mignani, Caraveo & Bignami 1997), 
no reliable spectral classification has been possible so 
far. The spectrum is dominated by a broad, asymmet- 
ric, Ha line (Duncan, Stewart & Haynes 1993; Mignani, 
Caraveo & Bignami 1997) with no prominent stellar emis- 
sion/absorption features (Johnston, Fender and Wu 1999). 
The radio counterpart of Cir X-l was identified soon af- 
ter the X-ray detection (Clark et al. 1975) and found to 
be variable with the same 16.6 day period (Haynes et al. 
1978). The radio picture features a point source which is 
embedded in a ~ 5 x 10 arcmin 2 nebula and is apparently 
the origin of elongated structures, extending on both ar- 
cmin (Stewart et al. 1993) and arcsec scales (Fender et al. 
1998). Such a morphology has been interpreted as a jet- 
powered synchrotron nebula, making Cir X-l the second 
galactic X-ray binary, besides SS433, which shows both 
radio jets and a radio synchrotron nebula. 
The shape of the radio nebula, protruding southward and 
pointing to the center of the nearby (~ 25 arcmin) radio 
SNR G321.9-0.3, is reminiscent of a residual wake left be- 
hind as the point source moves away from the SNR (see 
Fig.l, reproduced from Fig.l of Stewart et al. 1993). This 
suggested a tantalizing scenario where Cir X-l is a run- 
away binary system, with the neutron star member ejected 
by the supernova explosion responsible for SNR G321.9- 
0.3 (Stewart et al. 1993) between 2 x 10 4 and 10 5 yrs 



2 



R.P. Mignani et al.: Cir X-l 



ago (Clark, Caswell and Green 1975). Such a scenario is 
consistent, within the uncertainties, with the estimated 
distances of both Cir X-l and of the SNR. From H I ab- 
sorption measurements (Goss & Mebold 1977), the dis- 
tance to Cir X-l is anywhere between 6.7 and 8 kpc, while 
the distance to the SNR ranges from 5.5 kpc (Case & 
Battacharya 1998) to 6.5 kpc (Stewart et al. 1993). If the 
above scenario is correct, for a SNR age T (yrs), Cir X- 
1 should move Northward in the plane of the sky with a 
proper motion \i ~ 15 (T/lQ 5 yrs)~ 1 mas yr -1 . This value 
would imply a tranverse velocity with respect to the cen- 
ter of the SNR of v ~ 5 \i dk pc km s _1 , where fi is in units 
of mas yr and dkpc is 

the SNR distance in units of kpc. 
From the combination of all the possible values of both 
the SNR age and distance, we would thus expect Cir X-l 
to move Northward with a proper motion ranging between 
15 and 75 mas yr -1 , corresponding to a transverse veloc- 
ity between 400 and 2400 km s _1 . 

Thus, the Cir X-l angular displacement is a simple, yet 
powerful, tool to prove or disprove its association with 
SNR G321.9-0.3. Measuring a Northward angular dis- 
placement > 15 mas yr _1 would secure the association, 
while the lack of such a displacement, or a displacement 
in a different direction, would rule it out. In principle, 
both the optical and radio counterpart could be used to 
gauge the precise position of the source. However, the dra- 
matic fading of the radio point source (Stewart et al. 1991) 
hampers a measurement in the radio band and leaves the 
burden of the proof to astrometry of the optical counter- 
part. 

The sharpness of view of HST makes it the ideal instru- 
ment to settle such a proper motion issue. Thus, we have 
taken advantage of an HST/WFC archival image of the 
field as well as of recently acquired WFPC2 ones to carry 
out the proper motion measurement. The observations 
and data reduction are described in Section 2 while the 
results are presented and discussed in Sections 3 and 4. 

2. Observations. 

2.1. The data sets 

The first image of the Cir X-l field was obtained with 
the pre-repair HST on October 20 th 1992 as a target ac- 
quisition image for spectroscopy observations with the 
Faint Object Spectrograph (Mignani et al. 1997). The 
observation was performed with the Wide Field Camera 
(WFC) (O'.'l/pixcl) through the wide band filter 785LP 
"I-band filter" (A = 9703 A; AA = 1669.5 A), where 
the source is brighter (m785Lp — 18). The exposure time 
was 300 s. A new observation of Cir X-l was performed 
on May 19 th 2001 with the Wide Field and Planetary 
Camera 2 (WFPC2), as a snapshot programme. In order 
to achieve the highest positional accuracy, the target was 
positioned at the center of the Planetary Camera (PC) 
chip (O'.'045/pixcl). For consistency with the first obser- 
vation, the 785LP filter was used. Three observations of 
300 s each were obtained to allow for a better cosmic ray 




Fig. 1. ATCA radio map of the Cir X-l region, adapted 
from Fig.l of Stewart et al. (1993). Cir X-l is at the center 
of the compact radio nebula recognizable near the center 
of the field, with SNR G321. 9-0.3 due south of it. 

filtering. After the standard HST pipeline reduction (de- 
biassing, dark removal, flatfielding) the 2001 images have 
been combined with a median filter and stacked (Fig. 2). 
Residual cosmic ray hits were filtered out using specific 
IRAF/STSDAS routines. 

2.2. The data analysis 

The cleaned images are the starting point for our astro- 
metric analysis. Following a robust method successfully 
applied in several previous astrometric works (see e.g. 
Caraveo et al. 1996; Caraveo & Mignani 1999; De Luca 
et al. 2000,2001; Mignani et al. 2000,2001: Caraveo et al. 
2001), this relies on an accurate image superposition which 
is obtained by optimizing a linear coordinate transforma- 
tion between two reference grids determined by a set of 
common reference objects. 

In order to obtain an accurate position determina- 
tion, the reference objects have to be carefully selected 
i.e. they must be (i) not too close to the CCD edges, (ii) 
not extended, (iii) not saturated but not too faint. In or- 
der to avoid problems with the inter-chip astrometry (see 
Caraveo & Mignani 1999), the sclccton of the reference 
objects was limited to the chip where the target star was 
imaged, i.e. WFC #2 in the 1992 observation and PC in 
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Fig. 2. WFPC2 image of the Cir X-l field taken in 2001 
with the WFPC2 through the F785LP filter. Only the PC 
field is shown. North is up, East is left. The two ticks mark 
the position of the companion star of the X-ray source. 

the 2001 one. In the rather crowded field of Cir X-l, 46 
good reference sources, matching the above constraints, 
could be found. Their positions were computed with a 2- 
D gaussian fitting of their intensity profiles using centering 
boxes optimized in order to obtain stable values for both 
the centroids coordinates and their uncertainties (see De 
Luca et al. 2001 for a detailed discussion of both statistical 
and systematic errors involved in the procedure). The ref- 
erence objects were positioned with an accuracy between 
0.05 and 0.08 pixel (i.e. between 5 and 8 mas) in the WFC 
image and between 0.05 and 0.1 pixel (i.e. between 2.3 and 
4.5 mas) in the PC one, depending on the object's S/N. 
The position of the target star was computed in the same 
way, with an accuracy of about 0.07 pixel (7 mas) and 
of about 0.06 pixel (~3 mas) in the WFC and in the PC 
image, respectively. We note that the uncertainties on the 
centroids in both the WFC and the PC images are some- 
what higher than the ones measured in other fields (see 
e.g. De Luca et al. 2000). This can be ascribed, at least in 
part, to the peculiar PSF in the long-pass F785LP filter, 
characterized by broad and relatively bright wings. This 
is the same effect observed in filter F1042M (it is known 
as the "F1042M PSF anomaly" ) and it is due to the CCD 
detector becoming transparent at longer wavelengths, so 
that light is reflected and scattered by the back of the 
CCD producing a defocused halo (see WFPC2 Instrument 
Handbook for further details). In addition, for the 2001 
snapshot observations the analysis of the spacecraft atti- 
tude showed a jitter of the telescope about 2 times larger 
than average (w 5 mas rms compared to 3 mas rms) and 



distributed along a preferential direction. Moreover, we 
note that the 2001 observations were performed during 
the day/night transition, which occurred during the first 
exposure. Thermal effects, such as "breathing" , could thus 
be present, determining a slight defocussing of the camera. 
All these effects clearly influence the PSF, thus affecting 
the centering accuracy. 

The derived sets of coordinates were then corrected for 
known systematics, like the "34*™ row defect" of the PC 
(Anderson & King 1999), and for the effects of the distor- 
tions of the WFC and the PC. While for the WFC only the 
geometric transformation of Gilmozzi et al. (1992) is avail- 
able, for the WFPC2 we have used the most recent one 
provided by Casertano & Wiggs (2000), which accounts 
for the evolution of the distortion map of the WFPC2 due 
to the drift of the instrument in the Optical Telescope 
Assembly focal plane. 

The superposition of the frames was performed using the 
PC image as reference. First, the PC image was aligned 
in Right Ascension and Declination, according to the 
telescope roll angle. Then, the WFC reference grid was 
registered on the PC one by fitting a linear coordinate 
transformation, accounting for 2 independent scale fac- 
tors, 2 translation factors and a rotation angle. Following 
Caraveo et al. (2001), we checked for systematics due, 
e.g., to unknown displacements of the reference objects 
by repeating the frame superposition 46 times excluding, 
in turn, one of our reference stars. Four objects yielding 
relatively high deviations (between 0.5 and 1 PC pixel 
per coordinate) were thus identified and rejected. Different 
centering algorithms (e.g. separate gaussian fitting on the 
X and Y intensity distributions) and coordinate trans- 
formation routines were also tested. In all cases we ob- 
tained statistically indistinguishable results, showing that 
our procedure is correct and free of systematics. 
The final frame registration yielded a rms values on the 
residuals of the reference objects' coordinates of 0.33 PC 
pixel (15 mas) in Right Ascension and 0.25 PC pixel 
(11.4 mas) in Declination. Such rms values are a fac- 
tor of 4 higher than those obtained by De Luca et al. 
(2000a) for a WFC-to-PC superposition in the Vela pul- 
sar field. While the higher uncertainties in the centroids 
determination addressed above play a role, they can not 
account entirely for the larger residuals obtained in the 
image registration. This effect might be related to the 
likely dependance on the filter wavelength of the trans- 
formation coefficients used in the geometric correction. 
Unfortunately, since such dependance has been charac- 
terized for the WFPC2 (Trauger et al. 1995) but not for 
the WFC, the accuracy of the image registration might be 
affected by a non-optimized geometric correction. 

3. Results. 

We have plotted in Fig. 3 the relative positions of our tar- 
get together with the residuals on the reference stars posi- 
tions. Since the difference between the positions of Cir X-l 
is comparable with the residuals measured for the 42 refer- 



4 



R.P. Mignani et al.: Cir X-l 



0.4 - 



-0.5 



0.5 



Relative position (pixel) 



Fig. 3. Relative positions (crosses) of the Cir X-l com- 
panion star measured at the epoch of our two observation. 
The error bars for the 2001 position (defining the 0,0 in 
the plot) account only for the centering uncertainty, while 
for the 1992 position they account also for the frame regis- 
tration uncertainty. We overplotted as dashed arrows the 
residuals on the reference stars coordinates after the im- 
age superposition. Their rms value represents the frame 
registration uncertainty (see text). The displacement of 
Cir X-l is comparable with the residuals of the reference 
stars. 



ence stars, we conclude that the Cir X-l counterpart does 
not show any significant displacement. The position offsets 
amount to 0.24±0.33 PC pixel in RA and -0.15±0.25 PC 
pixel in Dec. By applying the nominal plate scale tran- 
formation (PC pixel size 0'.'0455), these values translate 
into an angular displacement in the plane of the sky of 
-11±15 mas in RA and -7±11 mas in Dec (neglecting the 
correction for the Sun's proper motion). Although such 
displacements are statistically null, it is worth comparing 
the derived proper motion with the predicted values. Over 
the time interval of 8.6 years spanned by our HST obser- 
vations, we obtain a total yearly displacement of 1.5±1.6 
mas yr _1 , setting the 3cr upper limit on the proper motion 
to w 5 mas yr . Such a value rules out even the minimum 
value of m 15 mas yr^ 1 required by the proposed associa- 
tion with G321.9-0.3 (see Sect. 1). Moreover, we note that 
also the derived constraints on the possible proper motion 
directions (Fig. 3) are incompatible with the relative po- 
sitions of Cir X-l and the SNR center, which are roughly 
aligned North-South (see Fig. 1). Thus, we conclude that, 
contrary to earlier speculations, Cir X-l and SNR G321.9- 
0.3 are not associated. 



4. Discussion. 

Both the lack of a measurable proper motion for Cir X- 
1 and its resulting non-association with SNR G321.9-0.3 
have interesting astrophysical implications. Firstly, it is 
important to note that our result frees the neutron star 
age from the tight constraints imposed by the proposed 
SNR association. Indeed, an age greater than 100 000 years 
would be more easily compatible with the estimated time 
scales of the accretion-driven decay of the neutron star 
magnetic field (see, e.g., Konar & Bhattacharya 1997, 1999 
and reference therein) down to the values expected for the 
occurance of type I X-ray bursts (Lewin & Joss 1983). 
As far as the dynamics of the binary system are concerned, 
the null proper motion measurement sets a 3tr upper limit 
of 24 dkpc km s _1 on its relative tranverse velocity, where 
dkpc is the source distance in kpc. Such a value would 
imply that Cir X-l is not moving with an anomalous ve- 
locity larger than, e.g., 400 km s -1 unless its spatial veloc- 
ity is essentially radial. A possible indication for an high 
radial velocity (« 430 km s^ 1 ) was indeed found in the 
spectroscopic observations of the Cir X-l optical counter- 
part (Johnston, Fender & Wu 1999), although it could not 
be unambiguously confirmed by the more recent work of 
Johnston et al. (2001). Together with our constraint on 
the tranverse velocity, an accurate spectroscopic measure- 
ment of the radial velocity would be critical to determine 
the dynamics of the binary system and to compute the 
kick velocity impressed by the SN explosion (Tauris et 
al. 1999). Moreover, the above constraints on the binary 
system dynamics may require a different scenario for the 
source interactions with the surrounding environment. In 
particular, the lack of a proper motion in the expected 
direction prompts for an alternative explanation for the 
morphology of the arc-shaped radio features observed in 
the synchrotron nebula, so far interpreted as jets trailing 
behind by a fast-moving source. 

Last, but not least, our result calls, once again, for more 
caution in associating compact objects with SNRs purely 
on the basis of their proximity in the plane of the sky and 
on morphological arguments. 
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